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ABSTRACT 



> 

^f\ I Giant pulses are short, intense outbursts of radio emission with a power-law intensity 

distribution that have been observed from the Crab Pulsar and PSRB1937+21. We 
have undertaken a systematic study of giant pulses from PSRB1937+21 using the 

0\ '. Arecibo telescope at 430, 1420, and 2380 MHz. At 430 MHz, interstellar scattering 

broadens giant pulses to durations of -^ 50 ^secs, but at higher frequencies the pulses 



are very short, typically lasting only -^ 1-2 //sees. At each frequency, giant pulses are 
emitted only in narrow (< 10 /is) windows of pulse phase located ~ 55-70 /isec after 



H I the main and interpulse peaks. Although some pulse-to-pulse jitter in arrival times is 

^ ' observed, the mean arrival phase appears stable; a timing analysis of the giant pulses 

yields precision competitive with the best average profile timing studies. We have 
measured the intensity distribution of the giant pulses, confirming a roughly power-law 
distribution with approximate index of —1.8, contributing > 0.1% to the total flux at 
each frequency. We also find that the intensity of giant pulses falls off with a slightly 
steeper power of frequency than the ordinary radio emission. 



Subject headings: pulsars-pulsars: individual (PSRB1937+21, B0531+21) 



1. Introduction 

Despite their remarkable frequency stability, most radio pulsars exhibit considerable 
pulse-to-pulse intensity fiuctuations. A histogram of the pulse intensity typically shows a roughly 
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exponential or Lorentzian shape (e.g., Manchester &: Taylor 197^ ), with a tail extending out to 



perhaps ten times the mean pulse strength. 

The Crab pulsar (PSR B0531+21), by contrast, exhibits frequent, very strong radio pulses 
extending to hundreds of times the mean pulse intensity (e.g., Lundgren et al. 1995| ). These 



so-called "giant pulses" exhibit a characteristic power-law intensity distribution — a phenomenon 
that for many years was observed in no other pulsar. Despite hints from early observations 
( Wolszczan, Cordes, &: Stinebring 1984| ), it came as a surprise to find qualitatively similar behavior 



from PSRB1937-I-21, the first (and still fastest known) millisecond pulsar ( [Sallmen fc Backer 
1995i iCognard et al. 199"^ ). 



The Crab pulsar and PSRB1937-I-21 could hardly be less similar in their properties. The 
Crab pulsar, born in A.D. 1054, is the youngest known pulsar, and B1937-I-21 is one of the oldest, 
with a characteristic spin-down age of ~ 2 x 10^ yrs. The inferred surface magnetic field strength 
of the Crab pulsar, B = 4x 10^^ G, is about 10^ times as high as that of B1937+21. Although 
the Crab is the fastest of the high field pulsars, with a period of P = 33 ms, it is twenty times 
slower than the 1.56 ms B1937-I-21. The pulsars' only identified common feature is their similarly 
strong magnetic field strength at the velocity-of-light cylinder, ~ 10^ G, which is higher than that 



for any other known pulsar (Cognard et al. 1996). Whether this is coincidence, or in some way 



responsible for the observed giant pulse behavior, is unknown. 

All previous studies of giant pulses from B1937-I-21 have been at or near 430 MHz with the 
exception of some limited, inconclusive data at 1.4 GHz (Wolszczan et al. 1984). We present below 
the first multifrequency study of giant pulses in B1937-I-21. In §|2|, we describe the observations 
and giant pulse search algorithms. A comparison of the arrival times of normal and giant pulse 
emission is discussed in §|3[ Following this, we discuss pulse morphology in ^ and the limits of 
accurate timing in §|5[ Using the arrival data, we present intensity distributions of the giant pulses 
and the contaminant noise in §|6[ In ^ we discuss the approximate spectrum of the largest giant 
pulses. Finally, in ^ we point to open questions and necessary future observations. 



2. Observations and Signal Processing 

All observations were made using the Princeton Mark IV instrument (Stairs et al. 1999) at 



the 305 m Arecibo telescope, between 1998 February 21 and 1999 August 1. Observations were 
made at three frequencies, as part of an ongoing timing study of this pulsar. The pulsar signal 
strength varied because of interstellar scintillation; observations made during times of strong 
signal were retained for the giant pulse analysis. The data used for this work included 30 minutes 
(~ 10^ pulses) at 430 MHz, 4 hours (<- lO"^ pulses) at 1420 MHz, and 26 minutes (~ 10*^ pulses) at 
2380 MHz. 

After completion of the observations, the data were coherently dedispersed in software by 
convolution with a complex chirp function ( Hankins fc Rickett 1975| ), to remove the progressive 



phase delays as a function of frequency caused by free electrons in the interstellar medium (ISM). 



The analysis pipeline has been described by Stairs et al. 1999. At the dispersion measure of 



PSRB1937+21(71), with a 10 MHz bandwidth, the dispersive smearing at 430 MHz, 1420 MHz, 
and 2380 MHz totals 74.16 ms, 2.06 ms, and 0.44 ms, respectively, before coherent dedispersion. 
After dedispersion, the signals from the two orthogonal polarizations were squared and cross 
multiplied to produce the four Stokes parameters. Because we were analyzing data taken primarily 
for other purposes, in most cases insufficient calibration data were available for high precision 
polarization calibration. We therefore concentrate on analysis of the total intensity data. 

After coherent dedispersion, the standard timing analysis pipeline was used to fold the data 
synchronously with the known topocentric period of PSRB1937+21, to produce average profiles. 
In a parallel analysis, the data were searched for strong individual pulses. Initial exploratory 
analysis confirmed that all strong pulses were confined to fairly narrow windows on the tails of 
the main pulse (MP) and interpulse (IP). At 430 MHz, giant pulses were therefore identified as 
in Cognard et al. (1996), by measuring the integrated flux density in 150 /xsec windows located 
on the tails of the MP and IP. At higher frequencies, where the giant pulses were much narrower 
(many lasting only a few 0.1-0.2 //sec bins) and appeared in a region of pulse phase that was 
significantly wider than the individual pulses, giant pulses were identified by searching for pairs of 
bins with combined energy greater than a threshold level. Because of the relative computational 
efficiency of this procedure, we did not limit our search to particular regions of pulse phase. But 
as at 430 MHz, giant pulses were only detected in the tails of the MP and IP. (Note that this 
appears to be in conflict with the results of Wolszczan et al. (1984).) 

Interstellar scintillation strongly modulated the apparent intensity of the pulsar signal from 
day to day. Scintillation should affect the normal and giant pulse emission identically. Therefore, 
in order to compare giant pulse intensities observed on different dates, we accounted for these 
variations by calibrating the average normal emission flux density for each run to the power-law 
model given by Foster, Fairhead, and Backer (1991). They found that the flux density as a 
function of frequency could be expressed as F[mJy]= (25.9 it 2.6) j/-2.60±o.05^ where u is in GHz. 
We use units of [Jy] for flux density and [Jy-/isec] for integrated flux density. 



3. Giant pulse distribution in time and pulsar phase 

An important difference between giant pulses in the Crab pulsar and in PSRB1937-I-21 
appears in the distribution of the pulses with respect to the star's rotational phase. Individual 
Crab giant pulses arrive at phases distributed throughout most of the emission envelope of 
the normal pulse profile. In contrast, the early observations of giant pulses from B1937-I-21 at 
430 MHz found that they arrived only in narrow regions on the tails of the two normal pulse 
components ( [Backer 1995| , Cognard et al. 1996 ). With our greater sensitivity, we confirm this 



result at 430 MHz and find the same behavior at 1420 and 2380 MHz. 



Figures |l|-g| show the average flux due to all giant pulses as a function of pulse phase along 
with the average normal emission. For each run, giant pulses were selected using an intensity 
threshold chosen to minimize noise contamination. As is clear from these figures, the giant pulses 
occur well after the normal emission phase; in particular the giant pulses do not produce the 
"notch" emission on the trailing edge of the MP as had been speculated by Cognard et al. (1996). 

To characterize the average properties of the giant pulse emission, we have fit the average 
giant pulse profile at each frequency with a Gaussian model at 2380 and 1420 MHz and, to account 
for interstellar scattering, with a model consisting of a Gaussian convolved with an exponential tail 
at 430 MHz. We find best-fitting Gaussian widths (FWHM) of 3 fis (4 fis) for the MP (IP) giant 
pulse profile at 2380 MHz and 4.3 /xs (4.1 /xs) for the MP (IP) giant pulse profile at 1420 MHz. 
At 430 MHz, we find that the MP and IP giant pulse profiles can be adequately described as 
Gaussians with widths 6.6 //s and 6.4 ^s, respectively, convolved with an exponential scattering 
tail with T = 28 fxs. This scattering timescale r is similar to that estimated from measurements 
of the normal profile, confirming that the structure observed in the low frequency average giant 
pulse profile is dominated by propagation effects. (Because our exponential scattering model is 



probably an oversimplification of the true effects of scattering on the signal (e.g., Sallmen et al. 



1999]) , our estimates of the intrinsic width of the mean giant pulse profile at 430 MHz should be 



considered an upper limit.) 

At high frequency, the windows in which giant pulse emission occur are much narrower than 
the average pulse emission windows. Indeed, they are remarkably narrow both absolutely and as a 
fraction of the pulsar period, each corresponding to less than one degree of rotational phase. They 
are, we believe, the sharpest stable features ever detected in a pulsar profile. It is this property 
that makes the mean giant pulse emission from PSRB1937-I-21 a potentially valuable fiducial 
point for high-precision timing observations, as discussed in §|5| below. 

We have also investigated the distribution of giant pulses over time, finding consistency with 
Poisson statistics (so neighboring giant pulses appear uncorrelated) . The distribution in time 
during one observation is illustrated in Fig. 0, which displays all pulses with integrated flux density 
> 30 Jy-fisecs. Also shown are the brightest pulses, with integrated flux density > 80 Jy-/isecs. 
The arrival time uncertainty for any given giant pulse is ^ 0.5 /isecs; this will be discussed further 
discussed in [Q We find no difference in the distributions of the more powerful and less powerful 
giant pulses, except for rate. The giant pulse distribution also appears very stable over this 
half-hour scan, with no apparent drifting or nulling periods, although the density drops off slightly 
at the end of the scan, especially in the IP, as interstellar scintillation reduces the mean pulsar 
flux relative to our threshold values. 



4. Individual Pulse Morphology 

The average giant pulse profiles constrain the properties of the giant pulse emission region. 
Also of interest are the properties of individual giant pulses, which constrain the giant pulse 
emission mechanism itself. To study individual giant pulses, which are narrower than the mean 
giant pulse envelope, we must account for scattering effects even at frequencies above 430 MHz. As 
discussed above, scattering by a thin turbulent screen degrades a signal by effectively convolving it 
with a one-sided exponential tail. (Although more sophisticated models of the ISM will have more 
complex effects on the observed signal, we find the simple thin-screen model adequately describes 
the vast majority of giant pulses we have observed.) We have fit all candidate giant pulses with 
a Gaussian, ^exp{— (t — to)^/(2o'^)}, convolved with an exponential tail, e~*''^/r, with four fit 
parameters, A, to, a, and r. Shown in Fig. ^ are the strongest giant pulses found at each frequency 
along with their fitted convolved Gaussians. 

After visually inspecting many giant pulses at all three frequencies, we see no evidence for 
intrinsic multiple-peaked emission, which contrasts with giant pulses observed from the Crab 
pulsar ( Sallmen et al. 19991) . Although receiver noise dominates the profiles, most pulses show a 



fast rise followed by an exponential decay, consistent, on average, with scattering. We have more 
quantitatively verified our model for the pulse morphology by cross correlating all of the giant 
pulses with a standard giant pulse, shifting by the lag which maximizes the cross correlation (to 
correct for the observed pulse-to-pulse jitter), then folding. For any choice of the standard giant 
pulse, this has always produced a single-peaked, short-rise-time, exponential-decaying profile, 
which, itself, is well fit by our model. Gaussian widths (FWHM) of 6, 0.2, and 0.2 /isecs and 
scattering timescales of r ~ 29, 0.2, and 0.2 /usees were determined at 430, 1420, and 2380 MHz, 
respectively, for both the MP and IP shifted, folded giant pulse profiles. As previously mentioned 
in §& our exponential scattering model is probably an oversimplification of the true effects of 
scattering, implying that our estimates should be considered upper limits. 

Table |I] lists the range of scattering timescales, r [/isec], found at each frequency. Overall, 
we find the following ranges: r ~ 13-40 /is (430 MHz), r < 1.1 ^s (1420 MHz), and r < 0.4 fis 
(2380 MHz). These timescales are consistent with turbulent scattering, which has a frequency 
dependence of j/-^ to -4.4 (^g_g_^ Manchester fc Taylor 1977| ). In addition, we have determined 



approximate upper limits to the fitted Gaussian widths (FWHM) at each frequency of ^ 7 ^s 
(430 MHz), ^ 0.5 fis (1420 MHz), and ^ 0.3 /xs (2380 MHz), which have not, however, been 
included in Table [l|. 

Scattering also affects the normal emission, primarily observations at lower frequencies. The 
best-fitting scattering parameters were determined independently in the two pulse components. In 
the MP, our model finds r ~ 30 /isec, and in the IP, r ~ 40 ^usec. Both are in good agreement with 
the scattering time estimated from the giant pulse emission. The IP is well fit by this model, but 
our fit to the MP underestimates the fiux in the tail, probably implying the unresolved presence 
at this frequency of the "notch" feature that is resolved on the trailing edge of the MP at higher 



frequencies. (This is not unexpected, since the feature increases in strength relative to the MP 
peak with decreasing frequency.) Note that scattering not only broadens but also delays the peak 
of the low-frequency pulsar signal by an amount that depends on the pulse shape; variability of the 
scattering strength therefore introduces significant timing errors at low frequency. For the 1420 
and 2380 MHz normal emission, scattering is much less severe; therefore, we expect no significant 
scattering delay to the apparent MP and IP peak arrival times. 



5. Timing 

High-precision timing measurements benefit from a sharp-edged timing signal. For this 
reason, the timing properties of the giant pulses are of considerable importance, especially since 
very high signal-to-noise ratios can be obtained for mean giant pulse profiles by using a signal 
thresholding technique to eliminate data when no giant pulse is present. 

High-precision alignment of pulsar profiles at different frequencies is not straightforward when 
the pulse shape is variable, since the choice of fiducial reference phase for alignment is arbitrary. 
We have aligned the three profiles in Figs. |l]-^ by the peaks of the normal emission profiles, in the 
case of 430 MHz after accounting for the delay introduced by scattering. As is evident, at each 
frequency the giant pulses occur at approximately the same phase relative to the normal emission 
peaks. The slight delay at 2380 MHz and possibly at 430 MHz with respect to the 1420 MHz giant 
pulse profile is most likely due to our somewhat ad hoc alignment procedure. 

We present timing characteristics for each observation in Table IT]. Displayed for each scan are 
the separation in phase angle of the IP peak following the MP peak (Normal), the separation of 
the IP giant pulses from the MP giant pulses (Giant), and the delay (in [/is]) of the giant pulses 
with respect to the MP and IP emission peaks (with scattering taken into account for the 430 MHz 
observations, as will be discussed below). We now discuss each timing column in more detail. 

The separation of 57-58 /xsecs between MP peak and average giant pulse (though only a 
lower limit of 49 /usecs at 430 MHz), as well as that for the IP of 65-66 //sees, is the same at all 
three frequencies. This yields tight constraints on the relative geometry of normal and giant pulse 
emission regions. The separation between the MP and IP giant pulses is 189.5° at all frequencies, 
slightly larger than the 187.6° separation the MP and IP normal emission peaks (though only a 
lower limit of 185.6° can be quoted at 430 MHz). The individual pulse arrival phases are Gaussian 
distributed, with widths of cr = 1.5-2.0 //s, in good agreement with the width found for the average 
giant pulse profile (displayed in Figs. Q and |3|). This pulse-to-pulse jitter is evident in Fig. ^, 
where we have plotted the fractional giant pulse arrival bin versus pulse number for a 1420 MHz 
observation. 

It is interesting to ask whether observations of the giant pulse emission from PSRB1937+21 
can be used to carry out higher precision timing studies of the pulsar than have been possible using 
the relatively broad normal emission profile (e.g., Kaspi, Taylor, &: Ryba 1994| ). Using normal 



pulse timing techniques, absolute precisions as small as 0.12 fis have been obtained for B1937+21 
at 1420 MHz ( |Stairs 1998 ). Just as typical long term timing studies depend on long-term stability 



of the normal emission profile, timing studies using the giant pulses will depend upon long-term 
stability of the giant pulse emission phase distribution. The consistency of our results with those 
of Cognard et al. are encouraging, but careful observations at high frequency over a period of years 
will be needed to test the ultimate power of giant pulse observations for high precision timing. 

Nevertheless, we have done preliminary estimates of the obtainable timing precision using the 
current data set. As noted above, a single giant pulse at 1420 MHz can be used to estimate the 
pulsar phase to ^ 1.5 fis. The key question for timing is whether this precision can be improved by 
averaging over multiple pulses. If the pulse-to-pulse phase variations are uncorrelated, we expect 
the timing precision to improve as yJNg, where Ng is the number of consecutive points averaged. 
In Fig. |6| , we show the r.m.s. scatter a within individual days as a function of Ng. We find that 
timing precision improves as expected to the limit of our data sets, Ng = 64, corresponding to a 
level of 100-300 ns in a ten-minute observation. Also indicated in the figure is the best timing 
precision achieved using standard pulsar timing techniques, ~ 120 ns. Stairs et al. (1999) suggest 
that the limiting factor in their timing analysis may be pulse shape variations caused by variable 
interstellar scattering. If this is correct, then giant pulse observations may ultimately improve on 
normal pulse observations, because the effects of scattering can be much more easily measured 
and removed from the data. Whether these high-precision timing results on short timescales will 
lead to better long term timing depends primarily on the long term stability of the giant pulse 
emission characteristics, which will be studied in future work. 



6. Intensity Distribution 

The primary distinguishing characteristic of the giant pulses observed from the Crab Pulsar 
and PSRB1937+21 is, of course, their intensity, and particularly their extended power-law 
intensity distribution. In Figs. ^^, we plot the cumulative intensity distribution at each observing 
frequency. At low intensities, the distributions are dominated by the chi-square statistics of the 
noise and/or normal emission, but above a certain threshold the pulse strength distribution is 
roughly power-law distributed. Given our limited statistics for any given run, this simple model is 
adequate to describe the observed data. 

In Fig. ^, we plot the cumulative distribution of the integrated flux densities in 150 ^sec 
windows after the MP and IP during a single 15-minute observation at 430 MHz (MJD 51364). 
Also plotted is the —1.8 power-law slope, which Cognard et al. (1996) found for the cumulative 
giant pulse distribution at this same frequency. The mean signal-to-noise ratio for the normal 
emission peak in this observation was about 0.14. Contamination from this emission causes a 
noticeable deviation at low flux levels from what would be seen with giant pulses and receiver noise 
alone, causing a steepening of the distribution towards low integrated flux densities. Although 
the distribution of normal pulse signal strengths is not known, we expect that removing normal 



pulses would produce better accord with a single power law distribution for the giant pulses. In 
Fig. ^, we have plotted the cumulative distribution of all giant pulses found at 1420 MHz for all 
~ 4 hours of data (using the threshold detection algorithm described above). Again, generally 
power-law behavior is observed, with a similar power-law exponent around —1.8. In Fig. ^, we 
plot the cumulative intensity distribution for a ~ 26-minute 2380 MHz run (MJD 51391). Again 
we have plot a —1.8 slope for comparison, which appears to fit the MP giant pulses and the most 
energetic IP giant pulses well, though more data are needed to strengthen this result. 

We have also calculated the fraction, R, of the total pulsar emission at each frequency that 
emerges in the form of giant pulses. We find the following ranges at each frequency for this 
fraction: 0.15% < R < 9% (430 MHz), 0.13% < R < 4% (1420 MHz), and 0.10% < i? < 1% 
(2380 MHz), where the lower limits were determined directly from Figs. |^-^ and the upper limits 
were determined from Figs. |^-^ by assuming a cumulative distribution with power-law slope of 
—1.8 over all intensities for both the MP and IP giant pulses. From the folded normal emission 
alone, however, we can rule out large values for R in these calculated ranges, implying that the 
single power-law model may not be valid at low intensities. 

Also consistent with Cognard et al. (1996), we find significantly stronger giant pulses following 
the MP than the IP. At a given frequency (e.g., 1/minute), the ratio of the strongest giant pulse 
associated with the MP to the strongest associated with the IP is very roughly the same as the 
ratio of the peak flux density in the MP to that in the IP. Despite the fact that giant pulses are 
separated in pulse phase from the normal emission, this suggests a relatively close association 
between the emission processes. 



7. Spectrum 

The short timescale of the observed giant pulses imply that they are a relatively broadband 
phenomenon, with bandwidth greater than their inverse width. The limited bandpass available to 
the Mark IV instrument prevents stronger statements about the spectra of individual giant pulses. 
However, the similar arrival distributions and roughly similar arrival rates at each observational 
frequency point to a broadband phenomenon. 

Nevertheless, our observations can be used to constrain the average spectral properties of the 
giant pulse emission. To avoid complications arising from the use of different effective thresholds 
at each frequency (because of different source strengths and different receiver noise properties), we 
estimate the giant pulse spectrum by using the most powerful individual pulses observed during 
a given time period at each frequency. We have, as usual, calibrated our observations to the 
spectrum for the average normal emission flux density at each frequency from Foster et al. (1991), 
as discussed in ^. 



Figure IC shows the intensities in [Jy-/isec] of the top eight MP and top eight IP giant pulses 



at each frequency over 15 minutes, corresponding to the entire MJD 51364 (430 MHz) run and 



to 15-minute chunks from runs on MJD 50893 (1420 MHz) and MJD 51391 (2380 MHz) We find 
a somewhat steeper slope of —3.1 for the giant pulse spectrum, compared to the —2.6 slope for 
the normal emission spectrum. Although the precise slope of the giant pulse emission spectrum 
depends on the assumed normal emission spectrum, the result that the giant pulse emission is 
steeper than the normal emission is robust. However, if the apparent narrowing of the giant pulse 
emission region at higher frequency reflects a narrowing of a sharp emission cone, the slightly 
steeper spectrum of the giant pulse emission might be understood as the geometric effect of the 
position of the line of sight through the outer part of the emission region 

Assuming the giant pulses from the Crab pulsar are powered by curvature radiation, Sallmen 
et al. (1999) have calculated the necessary number and number density of radiating electrons. We 
perform the same calculation here for our observed giant pulses from PSRB1937+21. For coherent 
curvature radiation, the power emitted by A'^ electrons with relativistic factor 7 = (1 — w^/c^)~^' ^ 
travelling along magnetic field lines with radius of curvature pc is 

-— ^(^)^ (1) 

From the peak of the largest 1420 MHz giant pulse in Fig. ||, we can calculate the maximum 
number of electrons needed in one bunch to produce the observed emission (at frequencies greater 
than 430 MHz). Assuming this pulse is broadband with spectral index, —3.1, PSRB1937+21 is 
at a distance of 3.6 kpc, and giant pulse beaming is determined by the beam width 9 ~ 7~^, we 
find a total power at the peak of 7"^ x 10^^ erg-s~^. An upper limit on the giant pulse Gaussian 
FWHM of < 0.5 /xsec (at high frequencies), gives 7 > 500, which implies a substantially smaller 
power requirement than the total spin-down energy loss rate of 2 x 10^^ erg-s~^. Setting Pcurv 
equal to the observed power yields 

iv=ioi9('^vY^^y (2) 

V500/ VIO^ cm/ ^ ' 

In order to preserve coherence at the highest observational frequency, these electrons must 
fit within a cube with volume ^ A^, where A = 12.6 cm at 2380 MHz, implying an electon 
density of rig = 5 x ^^^^ {-^)~^ {j^^^^) cm~^, a value 50 times greater than the Goldreich- Julian 



density ( Goldreich fc Julian 1969| ), riQ^j = QB/2iTec = 10^^ (r^) '^ '^^^ ^' which is the electron 



number density required to power the normal emission via curvature radiation. Sallmen et al. 
(1999) similarly find a giant pulse electron density 100 times greater (for pc = 10^ cm) than the 
Goldreich-Julian value for the Crab pulsar. 



8. Discussion 

The discovery of giant pulses from a millisecond pulsar was unexpected, and early hopes that 
identifying commonalities between PSR 61937+21 and the Crab pulsar might lead to a better 
understanding of the giant pulse emission mechanism have so far not been realized. As our study 
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has confirmed, the giant pulse emission from PSRB1937+21 differs fundamentahy from the Crab 
giant pulses, despite their common power-law behavior. The most intriguing characteristic of 
the high-frequency pulses form PSRB1937+21 is their very narrow widths and the very limited 
regions of pulse phase in which they occur. 

Despite the continued mystery about their origin, it appears likely that giant pulses from 
PSRB1937+21 may prove a valuable tool. As we have discussed, their narrow intrinsic width 
and large flux make them attractive fiducial reference points for timing studies of a pulsar that 
is already among the most precisely timed (Kaspi et al. 1994). Another intriguing possibility is 
to use the giant pulses as bright flashbulbs to study scattering in the ISM. Because the intrinsic 
pulses are very narrow, the pulse shape as observed at the Earth traces out the time delays 
introduced by multipath scattering. The combination of this information with VLBI studies of 
the scattering disk is a potentially powerful tool for studying the three dimensional distribution of 
scattering material. 

It is important, of course, to identify giant pulse emission from other pulsars. As we have 
noted, the integrated giant pulse emission from PSRB1937+21 is too weak to make noticeable 
features in the average pulse profiles, so careful, single pulse studies are required. Very fine time 
resolution is needed to avoid substantially smearing the high-frequency pulses from B1937+21 and 
reducing their signal-to-noise ratio, and it is insufficient to study only the windows of pulse phase 
where normal emission is found, as has sometimes been the case in past studies with coherent 
dedispersion instruments. Instruments that use hardware dedispersion followed by sampling (like 
the Princeton Mark III, Stinebring et al. 1992| ) must also preserve sufficient dynamic range in the 



analog-to-digital conversion to detect and characterize pulses that are far stronger than the typical 
pulsar emission. Another subtlety concerning observation of giant pulses are possible projection 
effects. If the PSRB1937+21 giant pulse emission region is roughly Gaussian, then as the pulsar 
rotates, the giant pulse emission traces out less than 1% of the entire sky, though the likelihood 
of detecting giant pulses from a known radio pulsar might be substantially enhanced over this 
estimate by correlations between the angular patterns of the normal and giant emission. Although 
searches for giant pulse emission from slow pulsars have been unsuccessful, only a very small 
fraction of millisecond pulsars have been studied sufficiently to detect or rule out giant pulses. 
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on the draft. In addition, S.E.T. thanks his earlier collaborators in this work, I. Cognard and J. 
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Fig. 1. — Folded giant pulse emission (solid) and normal emission (dashed) at 2380 MHz (top), 
1420 MHz (middle), and 430 MHz (bottom) with unit intensities of 0.07 Jy, 0.24 Jy, and 3.0 Jy 
(calculated using the power-law model in Foster et al. (1991)) and total integration times at each 
frequency of 26 minutes, 4 hours, and 30 minutes, respectively. One pulse period equals 4096 bins 
with each bin ~ 0.38 ^usecs. The 2380 MHz profile has been rotated so that the MP peak occurs in 
the same bin as for the 1420 MHz profile. The 430 MHz profiles have been rotated so that the IP 
peak arrival bin (with scattering taken into account) is coincident with that of the 1420 MHz IP 
peak. The giant pulse profiles have been amplified by the stated amounts. (Because of differences 
in receiver noise levels and the giant pulse selection threshold, a comparison of the amplification 
factors is not physically meaningful.) 
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Fig. 2. — Average giant pulse emission (solid) and normal emission (dashed) for the main-pulse 
region. All other details are the same as for Fig. |l|. 



0.4 

0.2 





S 0.4 



<u 0.2 





0.4 

0.2 





J ftjfesJ^LfK 



H — I — I — I — I — — h 



"1 ' ' ' ' I r 

xl60 



'■'.V^ i>;"""','"/. - -'-"-;'.i' fet>iiV 



3200 



H — I — I — h- ^ — I — I — I — h 
xlOO 




d — ^ — I — ^ — ^ — I — ^ — I — ^ — I — ^ — ^ — I — I — ^ 



J I , , , , I L 

3300 3400 
Bin 



x250 




3500 



Fig. 3. — Average giant pulse emission (solid) and normal emission (dashed) for the interpulse 
region. All other details are the same as for Fig. |l|. 
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Fig. 4. — Distribution of giant pulses at 1420 MHz (on MJD 50865) with integrated flux density 
> 30 Jy-/xsecs, in both pulse phase (1 bin ~ 0.38 /isecs) and pulse number. Arrival uncertainty 
for any given giant pulse is ^ 0.5 /xsec. The left panel shows the MP region, the right panel the 
IP region. In each region, the pulses are confined to a very narrow window — 10 bins (of 4096 
total) corresponds to about 0.9° of pulse phase — and have a temporal distribution consistent with 
a Poisson process. The same is true for a subset of the strongest giant pulses, with integrated flux 
density > 80 Jy-^secs, which are indicated with open squares. 
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Fig. 5. — Strongest individual giant pulses at 2380 (top), 1420 (middle), and 430 (bottom) MHz. 
Flux density (solid) and fitted convolved Gaussians (long-dash) are shown. Each bin corresponds 
to ~ 0.38 fisecs. The sampling interval for the top and bottom plots is 0.2/isecs, and for the middle 
plot it is 0.1 /isecs. The 430 MHz data has been boxcar smoothed by 8 data points. 
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Fig. 6. — Estimated giant pulse timing precision. At each epoch where 1420 MHz observations 
were available, the standard deviation a of the giant pulse arrival phase is calculated for MP giants 
(triangles) and IP giants (X's), without averaging [Ng = 1). Then consecutive points are pairwise 
averaged [Ng = 2) and a is calculated again. Pairwise averaging continues as long as the number 
of points remaining is greater than one. The solid line with slope —1/2 is expected for uncorrelated 

— 1/2 

Gaussian timing noise, i.e., a (x Ng . We also plot the number of detectable giant pulses expected 
in 1-minute and 1-hour intervals for the MP (solid tick) and IP (dashed tick) . The dotted line is the 
timing precision, a ~ 120 ns, achieved with the Mark IV instrument using standard average-profile 
techniques (Stairs et al. (1999)). 
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Fig. 7. — Cumulative intensity distribution for all 150 ^sec windows after the MP (triangles) and 
IP (X's) at 430 MHz (15 min.-run on MJD 51364), where / denotes the fraction of pulses containing 
emission > S* in the specified windows. Because the giant pulses arrive within the normal emission 
window, the giant pulses are augmented by the normal emission. This steepens the observed giant 
pulse power law, leading to a poor fit to the displayed —1.8 power law. 
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Fig. 8. — Cumulative intensity distribution for candidate giant pulses at 1420 MHz using the entire 
4 hour data set, where S is calculated over 2 /isecs for each candidate and / denotes the fraction of 
pulses containing candidates arriving in the specified 11.4-^sec MP (triangles) or IP (X's) windows 
with intensity > S. Also shown is a power law with index —1.8. 
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Fig. 9. — Cumulative intensity distribution for candidate giant pulses at 2380 MHz (26-niin. run 
on MJD 51391), where S is calculated over 2 //sees and / denotes the fraction of pulses containing 
candidates arriving in the specified 15.2-//sec MP (triangles) or IP (X's) windows with intensity 
> S. The dotted line denotes all other candidates from the rest of the pulsar phase scaled to the 
giant pulse window width of 15.2 /xsecs. Also shown is a power law with index —1.8. 
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Fig. 10. — Giant pulse integrated flux density versus frequency for the eight strongest giant pulses 
in the MP and in the IP taken from equal-length scans of 15 minutes on MJD 51364 (430 MHz), 
50893 (1420 MHz), and 51391 (2380 MHz). At each frequency, the absolute flux density cahbration 
was calculated independent of scintillation using the spectral model of Foster et al. (1991). This 
power-law model has a spectral index for normal emission of —2.6; that slope is indicated above. 
There is some evidence that the spectrum of the strongest giant pulses is steeper; the best flt 
power-law model, with index —3.1, is also shown. 
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Table 1. 

Freq. Min. Normal Giant MP delay [^scc] IP delay [/isec] r [/isec] 



2380 


25.8 


187.9 ±0.3° 


189.5 ±0.5° 


58 ±1 


65 ±2 


^0.4 


1420 


242.6 


187.58 ±0.01° 


189.5 ±0.2° 


57 ±1 


65 ±1 


^1-1 


430 


30.0 


> 185.6 ±0.5° 


189.4 ±0.6° 


>49±2 


66 ±2 
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